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Abstract 

The new model of nn transitions in the medium based on unitary ^-matrix is consid- 
ered. The time-dependence and corrections to the model are studied. The lower limit on 
the free-space nn oscillation time r is in the range 10^^ yr > r > 1.2 • 10^ s. 
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1 Introduction 



In the standard calculations of ab oscillations in the medium [1,2] the interaction of particles a 
and b with the matter is described by the potentials Ua,b (potential model). Imf/f, is responsible 
for loss of 6-particle intensity. We consider the nn transitions [3,4] in a medium followed by 
annihilation: 

n^n^ M, (1) 

where M are the annihilation mesons. For the process (1) potential model is used as well [5-12]. 

In [10,11] it was shown that one-particle (potential) model mentioned above does not de- 
scribe the process (1) and thus total neutron-antineutron transition probability: the process 
(1) probability is 1^ ~ F (see (15)), whereas the potential model gives W ~ 1/F, where F is 
the annihilation width of n in the medium (see eq. (6) of Ref. [5] or eq. (1.6) of Ref. [9] or 
eq. (16) of present paper]). In the potential model the effect of final state absorption (annihi- 
lation) acts in the opposite (wrong) direction, which tends to the additional suppression of the 
nn transition. Since the annihilation is the main effect which defines the speed of process (1), 
the potential model should be rejected. This is because the unitarity condition is used for the 
essentially non-unitary S'-matrix [10,11]. The interaction Hamiltonian contains the antineutron 
optical potential Un and \vaUn plays a crucial role. The iS-matrix should be unitary. 

More formally, the basic equation Z^/^^i I ^/j P~ '^ImTa, S = 1 + iT follows from the 
unitarity condition [SS^)fi = dfi. However in the potential model the 5'-matrix is essentially 
non-unitary {SS~^)fi = 6fi + aji, resulting in | Tfi 2ImTii + an ^ llmTn because 

2ImTii is extremely small: 2ImTii < 10~^^ [10,11]. The above-given basic equation is inappli- 
cable in this case. 

The potential model was developed in [5-8]. In more recent papers the verious details of the 
model have been refined, in particular the parameters of optical potential. We don't dwell on 
these papers since the heart of the problem is in the non-Hermiticity of the optical potential. 

In [13] we have proposed the model based on the diagram technique for direct reactions which 
does not contain the non-Hermitian operators. Subsequently, this calculation was repeated in 
[14,15]. However, in [16] it was shown that this model is unsuitable because the nn transition 
takes place in the propagator which is wrong. The neutron of the nucleus is in the bound state 
and so the neutron line entering into the nn transition vertex should be the wave function, 
but not the propagator, as in the model based on the diagram technique. For the problem 
under study this fact is crucial. It leads to the cardinal error for the process in nuclei. The nn 
transitions in the medium and vacuum are not reproduced at all. If the neutron binding energy 
goes to zero, the result diverges (see eqs. (18) and (19) of ref. [13] or eqs. (15) and (17) of ref. 
[14]). So we abandoned this model [16] (for more details, see [17]). 
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In [10] the model which is free of drawbacks given above has been proposed (model b in the 
notations of present paper). However, the consideration was schematic since our concern was 
only with the role of the final state absorption in principle. In sect. 2 this model as well as the 
model with bare propagator are studied in detail. The corrections to the models (sect. 3) and 
time-dependence (sect. 4) are considered as well. In sect. 5 we sum up the present state of the 
investigations of the nn transition problem. 

Also we want to attract the attention to the problem under consideration because of the 
following reasons: 1) For the lower hmit on the free-space nn oscillation time the range of 
uncertainty is too wide (see eq. (38)). 2) The problem involves a number of questions which 
are of independent interest. 

The basic material is given in sects. 2 and 5. Figure 1 and analysis of amplitudes and 
Mb are of major importance. 

2 Models 

First of all wc consider the antineutron annihilation in the medium. The annihilation amplitude 
Ma is defined as 

</0| rexp(-i J dxn(x)) - 1 |Onp>= N (27:^ 5^ (p f - pi) Ma. (2) 

Here is the Hamiltonian of the n-medium interaction, | Oup > is the state of the medium 
containing the n with the 4- momentum p = (e, p); < / | denotes the annihilation products, N 
includes the normalization factors of the wave functions. The antineutron annihilation width 
r is expressed through Ma'. 

jd^\Ma\^ . (3) 

For the Hamiltonian % we consider the model 

H{t) - j d'xUix) = Ha{t) + V, (4) 

where Ha is the effective annihilation Hamiltonian in the second quantization representation, 
V is the residual scalar field. The diagrams for the model (4) are shown in fig. 1. The first 
diagram corresponds to the first order in Ha and so on. 

Consider now the process (1). The qualitative process picture is as follows. The free-space 
nn transition comes from the exchange of Higgs bosons with the mass niH > 10^ GeV [4] and so 
the subprocess of nn conversion is scarcely affected by a medium effects. /,From the dynamical 
point of view this is a momentary process: Tc ~ l/nin < 10^^^ s. The antineutron annihilates 
in a time Ta ~ l/F. We deal with two-step process with the characteristic time Tch ~ Ta- 
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Figure 1: Antineutron annihilation in the medium. The annihilation is shown by a circle 
The neutron wave function is 

n{x) = Q^^^^ exp{—ipx). (5) 

Herep = (e, p) is the neutron 4- momentum; e = p^/2m+C/„, where C/„ is the neutron potential. 
The interaction Hamiltonian has the form 

Hi = Hnn + H, (6) 

Hnn(t) = J d^x(enn^n(x)^n(x) + H.C.) (7) 

Here Hnn is the Hamiltonian of nn conversion [6], enn is a small parameter with enn = l/^nn, 
where is the free-space nn oscillation time; m„ = nin = ni. In the lowest order in Hnn the 
amplitude of process (1) is uniquely determined by the Hamiltonian (6): 

M = ennGoMa, (8) 
e-p72m-C/„ + iO' 

Pn — p, Cfi — e. Here Gq is the antineutron propagator. The corresponding diagram is 
shown in fig. 2a. The annihilation amplitude Ma is given by (2), where H — Ha + V^n^n. 
Since Ma contains all the n-medium interactions followed by annihilation including antineutron 
rescattering in the initial state, the antineutron propagator Go is bare. Once the antineutron 

annihilation amplitude is defined by (2), the expression for the process amplitude (8) rigorously 
follows from (6). For the time being we do not go into the singularity Gq ~ 1/0. 

One can construct the model with the dressed propagator. We include the scalar field V in 
the antineutron Green function 

G. = G. + G„VG. + ...= (j^^ = -i = -i, (10) 
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Figure 2: a nn transition in the medium followed by annihilation. The antineutron annihilation 
is shown by a circle, b Same as a but the antineutron propagator is dressed (see text) 

T: — V, where E is the antineutron self-energy. The process amplitude is 

M = EnnGdM^,, (11) 

GdM^, — GoMa (see fig. 2b). The block in the square braces shown in fig.l corresponds to the 
vertex function M^. The models shown in figs. 2a and 2b we denote as the models a and b, 
respectively. 

In both models the antineutron propagators don't contain the annihilation loops since the 
annihilation is taken into account in the amplitudes and M^; the interaction Hamiltonians 
Hi and unperturbed Hamiltonians are the same. If E — > 0, the model b goes into model a. In 
this sense the model a is the limiting case of the model b. 

We consider the model b. For the process width F;, one obtains 

2 2 

Ni j d^\M |2= jd^\Mb |2= ^F', (12) 

r' ^ Ni J \Mb\'^, (13) 

where F' is the annihilation width of h calculated through the (and not Ma). The nor- 
mahzation multiplier A^i is the same for F^ and F'. The vertex function is unknown. (We 
recall the antineutron annihilation width F is expressed through the amplitude Ma.) For the 
estimation we put 

M, = Ma, r' = F. (14) 

This is an uncontrollable approximation. 

The time-dependence is determined by the exponential decay law: 

2 2 

Wb{t) = 1 - e-"""' ^ ^r't = ^Tt. (15) 
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Equation (15) illustrates the result sensitivity to the value of parameter S. 

On the other hand, for nn transitions in nuclear matter the potential model gives the inverse 
F-dependence [6-11] 

r/2 4f^-/ 

where Un is the antineutron optical potential. The wrong F-dependence is a direct consequence 
of the inapplicabihty of the model based on optical potential for the calculation of the total 
process probability [11]. (The above-mentioned model describes the probability of finding an 
antineutron only.) 

Comparing with (15), one obtains 

Wh F^ 

where the values F = 100 MeV and E = 10 MeV have been used. The parameter E is uncertain. 
We have put E = ReC^ — C/„ 10 MeV only for estimation. 

The model b leads to an increase of the nn transition probabihty. The lower hmit on the 
free-space nn oscillation time increases as well. Let and Tpot be the lower limits on the 
free-space nn oscillation time obtained by means of eqs. (15) and (16), respectively; T^n is the 
oscillation time of neutron bound in a nucleus. The relationships between r**, Tpot and Tnn are 

r 

/ = V^Tpot = ^T-pot, (18) 

r' = ^V^' (19) 



where the well-known equation 



Tpot = 2JT^n/r (20) 



has been used. For estimation we take F = 100 MeV, E = 10 MeV and Tpot = 2.36 ■ 10^ s [18]. 
(The limit Tpot = 2.36 • 10® s was derived [18] through eq. (20) from experimental bound on the 
neutron lifetime in oxygen Tnn > 1.77 • 10^^ yr obtained by Super-Kamiokande collaboration 
[18].) Equation (18) gives 

= 5Tpot = 1.2 • lO'^ s. (21) 

If E — >■ 0, eq. (18) diverges: — >■ oo. This circumstance should be clarified; otherwise the 
model under consideration can be rejected. As we will see later, the correct formulation of the 
problem (on the finite time interval) leads to finite result for the model a, which justifies our 
approach. 

Wc return to the model shown in fig. 2a. We use the basis (n, n). The results do not depend 
on the basis. A main part of existing calculations have been done in n — n representation. 
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The physics of the problem is in the Hamiltonian. The transition to the basis of stationary 
states is a formal step. It has a sense only in the case of the potential model H = Hp^t = 
KeUfi — Un — i^f^ =const., when the Hamiltonian of n- medium interaction is replaced by the 
effective mass H — > Hp^t — nies because the Hermitian Hamiltonian of interaction of the 
stationary states with the medium is unknown. Since we work beyond the potential model, the 
procedure of diagonalization of mass matrix is unrelated to our problem. 
The amphtude (8) diverges 

M^ennGoMa^^^. (22) 

(See also eq. (21) of ref. [16].) These are infrared singularities conditioned by zero momentum 
transfer in the nn transition vertex. (In the model b the effective momentum transfer qq = 
V = takes place.) 

For solving the problem the field-theoretical approach with finite time interval [19] is used. It 
is infrared free. If if = Hpot, the approach with finite time interval reproduces all the results on 
the particle oscillations (see sect. 5.2 of ref. [12]). This is the test of above-mentioned approach. 
However, our purpose is to describe the process (1) by means of Hermitian Hamiltonian because 
the potential model describes the absorption wrongly. 

For the model a the process (1) probability was found to be [12,16] 

Wait)^Wfit) = el^f, ri»l, (23) 

where Wf is the free-space nn transition probability. Owing to annihilation channel, Wa is 
practically equal to the free-space nn transition probability. If i — > oo, eq. (23) diverges just 
as the modulus (22) squared does. If E — > 0, eq. (15) diverges quadratically as well. 

The explanation of the i^-dependence is simple. The process shown in fig. 2a represents 
two consecutive subprocesses. The speed and probability of the whole process are defined by 
those of the slower subprocess. If l/F <^ t, the annihilation can be considered instantaneous. 
So, the probability of process (1) is defined by the speed of the nn transition: Wa ~ Wf ~ t"^. 

Distribution (23) leads to very strong restriction on the free-space nn oscillation time [12,16]: 

r" = 10^*^ yr. (24) 

If E — > 0, Wb rises quadratically. So and r" can be considered as the estimations from 
below and above, respectively. 

3 Corrections 

We show that for the nn transition in medium the corrections to the models and additional 
baryon-number-violating processes (see fig. 3) cannot essentially change the results. First of 



7 



all we consider the incoherent contribution of the diagrams 3. In fig. 3a a meson is radiated 
before the nfi transition. The interaction Hamiltonian has the form 



Hj = J + H^- + H. 



(25) 



In the following the background neutron potential is omitted. The neutron wave function is 
given by (5), were p — (po, p) and po — -01, + p^/2m. 
For the process amphtude M30 one obtains 



M^,^gGer,nGM^^~'-\ (26) 

Po - ?o - "T. - (p - q)2/2m + iO' 
where q is the 4-momentum of meson radiated, M^""^) is the amplitude of antineutron annihi- 
lation in the medium in the {n — 1) mesons. As with model a, the antineutron propagator G 
is bare; the n self-energy S = 0. (The same is true for figs. 3b-3d.) 



n n n 



n n, n n 







n 



n 







n 



m - 

A 











Figure 3: Corrections to the models (a and b) and additional baryon-number- violating processes 
(c and d) 

If g ^ 0, the amphtude M^a increases since G ^ Gg, 

Gs = -. (28) 

Po — m — p'^/2m 

(The limiting transition g ^ for the diagram 3a is an imaginary procedure because in the 

vertex n — ?> n$ the real meson is escaped and so go > m,^.) The fact that the amplitude 
increases is essential for us because for fig. 2a g = 0. Due to this Gq ~ 1/0 and Wa ^ Wb. 
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Let and T^"'^ be the widths corresponding to the fig. 3a and annihilation width of n in 
the (n) mesons, respectively; T = Y.{n) Taking into account that F^"^ is a smooth function 
of y/s and summing over (n), it is easy to get the estimation: 

2 2 

Tsa fti 5 • 10-V%r ^ %r. (29) 
The time-dependence is determined by the exponential decay law: 

Wsa(t) ^ Tsat = (30) 

Comparing with (15) we have: Wsa/Wf, — V'^ /m% <^ 1. So for the model b the contribution of 
diagram 3a is negligible. 

For the model a the contribution of diagram 3a is inessential as well. Indeed, using eqs. 
(30) and (23) we get 

wut) r 



Wait) mir 

where we have put m$ = mTr. Consequently, if 



(31) 



mlt/r » 1, (32) 

and 

rt > 1 (33) 

(see (23)) then the contribution of diagram 3a is negligible. For the nn transition in nuclei these 
conditions are fulfilled since in this case F ~ 100 MeV and t = Tq = 1.3 yr, where Tq is the 
observation time in proton-decay type experiment [18]. In fact, it is suffice to hold condition 
(33) only because it is more strong. 

In the calculations made above the frcc-spacc nn transition operator has been used. This 
is impulse approximation which is employed for nuclear /3 decay, for instance. The simplest 
medium correction to the vertex (or off-diagonal mass, or transition mass) is shown in fig. 3b. 
In this event the replacement should be made: 

Cnn Cm = e„n(l Ae), (34) 

Ae = esb/^nn, where 635 is the correction to enn produced by the diagram 3b. For the model a 
the limit becomes 

= (1 + Ae)10^Vr- (35) 

Obviously, the Ae cannot change the order of magnitude of since the n — )■ n operator is 
essentially zero-range one. The free-space nn transition comes from the exchange of Higgs 
bosons with the mass niH > 10^ GeV [4]. Since nin ^ mw {jnw is the mass of W^-boson), the 
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renormalization effects should not exceed those characteristic of nuclear (3 decay which is less 
than 0.25 [20]. So the medium corrections to the vertex are inessential for us. The same is true 
for the model b. 

Consider now the baryon-number-violating decay n — > n$ [21] shown in fig. 3c. It leads to 
the same final state, as the processes depicted in figs. 2, 3a and 3b. Denoting | q [= g', for the 
decay width one obtains 



Qq = + The parameter e$ corresponding to the vertex n n$ is unknown and so no 
detailed calculation is possible. 

The baryon-number-violating conversion n — )■ A in the medium [21] shown in fig. 3d cannot 
produce interference, since it contains i^-meson in the final state. For the rest of the diagrams 
the significant interferences are unlikely because the final states in fiN annihilation are very 
complicated configurations and persistent phase relations between different amplitudes cannot 
be expected. This qualitative picture is confirmed by our calculations [22] for p-nuclear an- 
nihilation. It is easy to verify the following statement: if the incoherent contribution of the 
diagrams 3a-3c to the total nuclear annihilation width is taken into account, the lower limit on 
the free-space nn oscillation time r becomes even better. 

To summarise, the contribution of diagrams 3 is inessential for us. 

4 Time-dependence 

The non-trivial circumstance is the quadratic time-dependence in the model a: Wa ~ t"^- The 

heart of the problem is as follows. The processes depicted by the diagrams 2b and 3 are 
described by the exponential decay law. In the first vertex of these diagrams the momentum 
transfer (figs. 3a-3c), or effective momentum transfer (figs. 2b, 3d) takes place. The diagram 
2a contains the infrared divergence conditioned by zero momentum transfer in the nn transition 
vertex. This is unremovable perculiarity. This means that the standard S'-matrix approach is 
inapplicable [12,16,19]. In such an event, the other surprises can be expected as well. From this 
standpoint a non-exponential behaviour comes as no surprise to us. It seems natural that for 
non-singular and singular diagrams the functional structure of the results is different, including 
the time-dependence. The opposite situation would be strange. 

The fact that for the processes with g = the S'-matrix problem formulation (oo, — oo) is 
physically incorrect can be seen even from the limiting case H — 0: if Hi — Hnn (see (6)), the 
solution is periodic. It is obtained by means of non-stationary equations of motion and not 
jS-matrix theory. To reproduce the limiting case — > 0, i.e. the periodic solution, we have to 




(36) 
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use the approach with finite time interval. 

If the problem is formulated on the interval {t, 0), the decay width F cannot be introduced 
since F = J2f^i I >S'/j(oo, — oo) p /Tq, Tq — )■ oo. This means that the standard calculation 
scheme should be completely revised. (We would like to emphasise this fact.) The direct 
calculation by means of evolution operator gives the distribution (23). 

Formally, the different time-dependence is due to ^'-dependence of amplitudes. We consider 
eq. (26), for instance. If q decreases, the amplitude M^a increase; in the limit g — > it is 
singular (see (28)). The point q — corresponds to reahstic process shown in fig. 2a. The 
f^-dependence of this process is the consequence of the zero momentum transfer. 

The more physical explanation of the i^-dependence is as follows. In the Hamiltonian (25) 
corresponding to fig. 3a we put H — Hnn — 0. Then the virtual decay n — > n$ takes place. The 
first vertex of the diagram 3a dictates the exponential decay law of the overall process shown 
in fig. 3a. Similarly, in the Hamiltonian (6) corresponding to fig. 2a, we put H — 0. Then the 
free-space nn transition takes place which is quadratic in time: Wf{t) = e^^t^. The first vertex 
determines the time-dependence of the whole process at least for small F. We also recall that 
even for proton decay the possibihty of non-exponential behaviour is reahstic [23-25]. 

5 Discussian and summary 

In both models the antineutron propagators don't contain the annihilation loops since the anni- 
hilation is taken into account in the amplitudes Ma and Mf,. The alternative model containing 
the full in- medium propagator has been considered in [26]. The sole physical distinction be- 
tween models a and b is the definition of antineutron annihilation amplitude; or, similarly, the 
non-zero antineutron self-energy in the model b which is conditioned by residual scalar field. 
However, it leads to the fundamentally different results. 

If E — )■ 0, Wb{t) diverges quadratically. This circumstance should be clarified; otherwise 
the model b can be rejected. The calculation in the framework of the model a gives the finite 
result, which justifies our approach from a conceptual point of view and consideration of the 
model a at least as the limiting case. In fact, the model a seems quite realistic in itself. Indeed, 
we list the main drawbacks of the model b. 

1) The approximation — Ma is an uncontrollable one. The value of E is uncertain. 
These points are closely related. 

2) The diagram 2b means that the annihilation is turned on upon forming of the self-energy 
part — V (after multiple rescattering of n). This is counter- intuitive since at low energies 
[27,28] 

(Ja > 2.5(7,, (37) 
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where cr^ and ag are the cross sections of free-space nN annihilation and nN scattering, re- 
spectively. The inverse picture is in order: in the first stage of the n-medium interaction the 
annihilation occurs. This is obvious for the nn transitions in the gas. The model a reproduces 
the competition between scattering and annihilation in the intermediate state [26] . 

3) The time-dependence is a more important characteristic of any process. It is common 
knowledge that the t-dependence of the process probability in the vacuum and medium is 
identical (for example, exponential decay law (15)). In the model a the ^-dependencies in the 
vacuum and medium coincide: Wa ~ t^ and Wf ~ t^. The model b gives Wi, ~ t, whereas 
Wf ~ t'^. There is no reason known why we have such a fundamental change. 

H — Un, the model a reproduces all the well-known results on particle oscillations in 
contrast to the model b (see sect. 5.2 of ref. [12]). In other words, the model a reproduces all 
the results of potential model. (Recall that our purpose is to describe the process (1) by means 
of Hermitian Hamiltonian because the potential model describes the absorption wrongly.) 

The model a is free of drawbacks given above. The physics of the model is absolutely 
standard. For instance, for the processes shown in fig. 3 the antineutron propagators are bare 
as well. The same is also true for the diagram technique for direct reactions [13-15]. Besides, 
the process amplitude (8) is derived directly from the basic eqs. (6) and (2) which arc obvious. 

However, there is fundamental problem in the model a: the singularity of the amplitude 
(22). The approach with finite time interval gives the finite result, which justifies the models 
a and b at least in principle. Nevertheless, the time-dependence Wa ~ t^ and limit (24) seem 
very unusual. The corresponding calculation contains too many new elements. Due to this 
we view the results of the model a with certain caution. Besides, due to the zero momentum 
transfer in the nn-transition vertex, the model is extremely sensitive to the S. The process 
under study is unstable. The small change of antineutron self-energy E = — )■ E = 7^ 0, or, 
similarly, effctive momentum transfer in the nn transition vertex converts the model a to the 
model b: Wa — )■ with Wi, <^ Wa- (This is because the process amplitude is in the peculiar 
point (see (22)) owing to zero momentum transfer. For the processes with non-zero momentum 
transfer the result is httle affected by small change oi q.) Although we don't see the specific 
reasons for above-mentioned scenario, it must not be ruled out. This is a point of great nicety. 

Finally, the values = 1.2 • 10^ s and r" = 10^^ yr are interpreted as the estimations from 
below (conservative limit) and from above, respectively. So the realistic limit r can be in the 
range 

10^^ yr > T > 1.2 • lO'^ s. (38) 

Recall that for the free-space ab oscillations the ab transition probability is extremely sensitive 
to the difference of masses nia — nif, as well. For the neutron in the bound state {nn transition 
in finite nucleus) the result is the same [29], what is consistent with the qualitative process 
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picture: two-step process (dynamic nfi conversion [4], annihilation) with the characteristic time 
Tch ~ 1/r. The fact that process amphtude is in the pecuhar point (unremovable perculiarity) 
is the basic reason why the range (38) is very wide. The estimation from below = 1.2 ■ 10^ s 
exceeds the restriction given by the Grenoble reactor experiment [30] by a factor of 14 and the 
lower limit given by potential model by a factor of 5. At the same time the range of uncertainty 
of T is too wide. Further theoretical and experimental investigations are desirable. 
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